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Abstract. The galaxy cluster Abell 2218, at z = 0.171, is well-known for the discrepancy between mass estimates derived 
from X-ray and strong lensing analyses. With the present XMM-Newton observation, we are able to trace the gas density and 
temperature profiles out to a radius of ~ 1400 kpc (approximately the virial radius of the cluster). The overall surface 
brightness profile is well fitted over three orders of magnitude with a simple /3-model with a core radius of 0^95 and [3 = 0.63. 
The projected temperature profile declines steeply with radius (by ~ 50%), and is well described by a polytrope with parameters 
t = 8.09 keV and y = 1.15. The temperature map shows a pronounced peak in the central arcminute, where the temperature 
rises by a factor of two (from ~ 5 to ~ 10 keV). The mass profile, calculated assuming hydrostatic equilibrium and spherical 
symmetry, is best fitted with a King approximation to an isothermal sphere, implying a dark matter distribution with a central 
core, in contrast with the cusped cores found in more obviously relaxed clusters. The X-ray mass is approximately two times less 
than the strong lensing mass at r ~ 80 kpc, although the agreement between X-ray and weak lensing mass measurements 
at larger radius (r ~ 400 kpc) is slightly better. While the X-ray total mass estimates can vary by 30 per cent depending on 
the mass model, all measurements are significantly lower than the corresponding total mass from optical measurements. Given 
the X-ray results indicating considerable disturbance of the intracluster gas, leading to a probable violation of the assumption 
of hydrostatic equilibrium, and the observed substructure in the optical, suggesting a line-of-sight merger, it is unlikely that the 
different mass estimates of this cluster can be reconciled, at least with standard modelling assumptions. 
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1. Introduction 

In what has effectively become the standard scenario for the 
formation of structure in the Universe, galaxy clusters are born 
in the collapse of initial density fluctuations and grow under the 
influence of gravity. They are the nodes of cosmic structure, ac- 
cumulating mass through a combination of continuous accre- 
tion of smaller subunits along the filaments between clusters, 
interrupted by the occasional violent merger. Cosmological 
simulations show that, in a given cosmological volume, the 
sample of clusters will span the range from those which are 
relatively relaxed t o those which are clea rly in the midst of a 
major merger (e.g. jAscasibar et alll2003l) . not surprising given 
the hierarchical scenario. 

Given a sufficiently large sample of galaxy clusters, con- 
straints can be put on cosmological parameters. For instance, 
constraints on <x 8 and Q can be foun d by observing the 
cluster mass function and its evolution jPerrenodlll980i see 
iHenrvl I2004I for a recent review). Cluster masses can be 
estimated using traditional baryonic mass tr acers such as 
the v elocity dispersion of the galaxies (e.g., iGirardi et all 
Il998l) . or the X-ray emission of the intracluster medium (e.g., 
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iReiprich & Bohrineetl l2002l) . provided some assumptions are 
made about the relationship between the distribution of dark 
and luminous matter. For example, the fact that a given clus- 
ter is roughly spherical and in approximate hydrostatic equi- 
librium is a fundamental assumption when the total mass is 
derived from the X-ray data. On the other hand, gravitational 
lensing offers another view of the distribution of matter, allow- 
ing the calculation of the projected mass within a cluster, either 
from strong lensing at small scales or weak lensing on larger 
scales. However, because lensing measures the projected mass, 
strong lensing in particu lar is very sensitive t o additional matte r 
along the line of sight (Wambsgan ss. Bode & Ostr ikei| 12004 
and to the detailed dynamical state of the cluster llTorri et ; 
l2004l) . Indeed, it has been suggested that strong lensing arcs 
occur predominantly in cl usters which are dynamically m ore 
active than typical clusters jBartlemann & Steinmetd. [l996). 

Cluster mergers release vast amounts of energy into the 
X-ray emitting gas of the intracluster medium (ICM), leading 
to significant departures from hydrostatic equilibrium. Strong 
variations are produced in observable characteristics such as 
density, and, since much of the energy is dissipated in shocks, 
temperature and entropy. Temperature variations following a 
merger can survive ~ 4 - 6 times longer than density substruc- 
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Fig. 1. The left panel shows a large-scale XMM/DSS overlay image. The (log-spaced) contours are from the adaptively-smoothed, non- 
background subtracted, EMOS1+EMOS2 [0.5-2.0] keV image. The right panel shows an XMM/HST overlay image with the same contours. 



ture, and so the temperature distribution is a strong diagnostic 
of the dynamical state of a given cluster (e.g. jMarkevitch et all 
1998). 

A2218 was one of the first targets for a detailed mass study 
comparing X-ray and weak lensing data ( Sauir eset all ll~996 ). 
which concluded that these mass estimates were in agreement 
(see also lAllenl ll998). However, the strong lensing mass (es- 
timated using the remarkable arcs at a distance of ~ 80 h~^ 
kpc) was not in agreement wi th the X-ray estimate (e.g., 
iMiralda-Escud e & Babul 1 19951) . ROSAT studies suggested 
complicated X-ray stru c ture n ear the core dMarkevitc h, 1997; 
iNeumann & Bohringerl Il999l) : recent Chan dra observations 
have revealed this structure in greater detail ( Machacek et al., 
120021) . This has given rise to the general conclusion that 
the cluster is dynamically active, which is also supported 
by the presence of mass substructu re in subseq uent lens- 
ing studies (the latest being Smi th et alJ 12004), and in- 
deed substructure in the galaxy distribution in the optical 
(Girardi & Mezzetti, 2001). The recent Chandra temperature 
map o f the central regions s hows quite strong temperature vari- 
ations (Govoni et al., 2004). The cluster also hosts a weak radio 
halo (Giovannini & Ferretti, 2000). 

With the present XMM-Newton observations we derive gas 
density and temperature profiles out to large radius, generate a 
large scale temperature map, and, assuming hydrostatic equi- 
librium and spherical symmetry, calculate a mass profile of the 
cluster. As we will see, the assumptions inherent in the mass 
profile calculation are probably violated, at least in the cluster 
centre. 

The cluster lies at a redshift of z - 0.171. In our chosen 
cosmology (H = 70 km s -1 Mpc -1 , Q„, = 0.3 and Q A = 0.7), 
1' corresponds to ~ 175 kpc. Errors are quoted at the la level 
unless otherwise stated. 



2. Observations and data reduction 

A2218 was observed for ~ 18 ks on 2002 September 28. The 
observation showed no soft proton flares in the high energy 
(> 10 keV) light curves. There is a slight (~ 15 per cent) in- 
crease in the broad-band [0.3 -10.0 keV] count rate at the be- 
ginning of the observation, lasting ~ 2 ks. Since the increase 
is very weak and of short duration, and has no effect on the 
derived normalisation at high energy, we have used the entire 
observation. 

Except where noted, in the following analysis we use 
EMOS PATTERNS 1-12 and EPN PATTERN events. To 
correct for vignetting, the photon weighting method as im- 
plemented in the SAS task eviqweiq ht was applied to each 
event file. The blank-sky event lists of lRead & Ponmanl J2003) 
were used as background files, with these event lists being sub- 
jected to the same screening criteria, PATTERN selection and 
vignetting correction as the observation event files. 



3. Gas density distribution 

3.1. Morphology 

An XMM/DSS overlay image is shown in Fig.^ The cluster 
exhibits slightly elliptical, symmetric X-ray emission at large 
scale. Near the core, the X-ray emission becomes elongated in 
the SE-NW direction, which is the same orientation as the line 
joining the two brightest cluster galaxies. As can be seen in the 
XMM/HST overlay image, the peak of the [0.5-2.0] keV band 
X-ray brightness is slightly offset from the cD galaxy. Both 
the cD and the [0.5-2.0] keV X-ray peak are offset from the 
large-scale centre of symmetry. Since the complicated X-ray 
morphology of the central re gions has been exte nsively exam- 
ined using both ROSAT HRI jMarkevitchl 1 1 997l) and Chandra 
iMachacek et alJ il2002h . we will not discuss this further. 
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Table 1. /J-model fits for various sectors of A2218, errors are 90 per 
cent confidence. Sectors are defined anticlockwise from the positive 
Y-axis. 



Sector 


P 


r c 

(arcmin) 


* 2 /dof 


Full 


63 +0 01 


o:95 + °;° 2 

-0:02 


208.8/168 


0-90 

90-180 

180-270 


63 +0 01 

-0.01 

061 +o.oi 
74 +0 03 

'^-0.02 


i:o4+ o: 05 

-0,05 
, 98+ 04 

-o;o4 

i:21+ 07 
-0,06 


140.0/127 
275.6/148 
86.9/96 


270-360 


66 +0 01 


0:87 +0 ; 04 

-0.04 


161.8/115 



3.2. Surface brightness profile 

After identification and exclusion of 28 serendipitous sources, 
an azimuthally averaged surface brightness profile was gener- 
ated from source and background event files of each camera 
in the 0.3-3.0 keV band. Events were extracted in circular an- 
nuli centred on the emission peak of the cluster. Background 
subt raction for each c amera was undertaken as described in 
IPratt & Arnaudl J2002h . The EMOS and EPN profiles are in 
excellent accord and thus the individual camera profiles were 
coadded, and the total profile was then binned to a S /N ratio 
of 3cr above background. This profile is shown in Figure |2 
Cluster emission is detected out to slightly more than 11', cor- 
responding to slightly less than 2 Mpc for ou r chosen cos- 
molog y. For comparison, the ROSAT profile of ICannon et alJ 
( 1999) was detected out to ~ 12'. 

3.3. Density profile modelling 

In order to calculate the gas density profile, we fitted the sur- 
face brightness profile with analytic models. The fitting process 
uses simplex minimisation and takes into account projection 
effects in the line of sight. The models we use are convolved 
with the XMM-Newton PSF and binned into the same bins as 
the observed profile. A standard isothermal yS-model is an ac- 
ceptable fit to these data (x = 208.8 for 168 degrees of free- 
dom), and the values found for /? and r r are in good a greement 
with ROSAT-derived values (e.g., Birkinshaw & Hushes 1994; 
ISauires et al1ll996l iMarkevitchll 19971 ICannon et al.lll999h . In 
Figure |2 the best-fitting model is shown compared to the sur- 
face brightness profile. In Tablef^we give the /?-model param- 
eters for this fit as well as for surface brightness profiles ex- 
tracted in 90° sectors defined anticlockwise from the positive 
Y-axis. In all instances bar the 90° -180° sector, ay8-model is an 
adequate fit to the surface brightness. As a further check we cal- 
culated surface brightness profiles centred on t he broad-b and 
X-ray peak found in the Chandra analysis of Mac hacek et alJ 
(2002). These profiles yield best-fitting /?-model parameters 
which are in agreement with those in Tabled 

On close examination of Fig. |2] it can be seen that the /3- 
model is an excellent description of the data over three orders 
of magnitude in surface brightness. Related to this is the lack of 
evidence fo r a turnover in the sur face brightness at very large 
radius (c.f., Vikhli nin et al.lll999i) . although this latter point is 
highly sensitive to the exact background subtraction method. 
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Fig. 2. The combined EPIC 0.3-3.0 keV surface brightness profile 
of A2218. The profile is background subtracted and corrected for vi- 
gnetting. The solid line is the best fitting /?-model. See Table Q for 
model details. 

4. Spatially resolved spectroscopy 

In the following, the background spectra were normalised us- 
ing the count rates in the source-free [10.-12.] keV (EMOS) 
and [12.-14.] keV (EPN) bands. The normalisation factors used 
were 1.09, 1.11 and 1.15 for EMOS1, EMOS2 and EPN, re- 
spectively. 

4.1. Global spectrum 

We extracted a global spectrum from all events within a circle 
of radius 5.'1 to allow for direct comparison with the Chandra 
results. We fitted the data using a simple single temperature 
absorbed MEKAL model with the redshift fixed at the average 
redshift of the cluster (z — 0.17). 

It was found that naive fitting of the spectra in the band 
[0.3-10.0] keV with the absorption fixed at the galactic level 
(3.24 x 10 2() ctrT 2 ) produced results where the EPN temper- 
ature was significantly lower than that found by the EMOS. 
Conversely, if the absorption was let free, then the EPN ab- 
sorption was significantly lower than that found by the EMOS 
(which was in good agreement with the galactic value). We 
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thus a dopted the ad hoc strategy discussed in IPratt & Arnaudl 
fixing the low-energy cutoff at 0.6 keV for the EMOS 
and 1 .0 ke V for the EPN, which leads to good agreement be- 
tween EPIC temperatures with the absorption fixed at the galac- 
tic value. 

The global temperature in this case is 6.63 + 0.27, with 
an abundance of 0.13 + .04 solar (90% confidenc e), the 
latter measured relative to lAnders & Grevessel Jl989l) . This 
is in good agreement with pre viously determ i ned va lues, 
such as the Ginga analysis of iMcHardv et all dl990|). the 
ASCA an alyses of iMushotzky & Loewensteinl Jl997l)Kller] 
i 19981) andlCannon et alJ i 19991) . and the Chandra analysis of 
Machacek et al.N2002l). 



4.2. Radial temperature profile 

We then calculated a simple projected radial temperature pro- 
file. After exclusion of point sources, spectra were extracted in 
10 annuli centred on the X-ray emission peak and extending to 
the radius of maximum detection of the surface brightness pro- 
file (i.e., 0'- 11'). The spectra were binned to 3<x to allow the 
use of Gaussian statistics. EPIC spectra were then simultane- 
ously fitted in the [1.-10.] keV band with a single temperature 
MEKAL model with free abundances and the absorption fixed 
at the galactic value. We choose a relatively high low-energy 
cutoff to avoid problems related to the variation of the soft X- 
ray background with position on the sky. 

While we have already normalised the background using 
the count rates in the source-free high energy bands, back- 
ground subtraction is critical in spectroscopic studies at the 
lowest surface brightness levels, and can have significant ef- 
fects on the radial profile at high radii. We thus re-fitted the 
spectra with the background at +10 per cent of nominal, the 
variation seen in the typical quiescent background level. The 
temperatures of all annuli agree within the lcr error bars what- 
ever the normalisation. 

The projected temperature profile is shown in Fig. [3] The 
error bars are the upper and lower limit temperatures for each 
annulus taking into account variations of +10 percent in the 
background normalisation. These uncertainties are the total 
systematic plus statistical errors and should be viewed as very 
conservative estimates. As expected, the external annuli ex- 
hibit the largest uncertainties; however, the robustly determined 
lower limit of the temperature in the final annulus is 2.3 keV. 
The temperature profile is in excell ent agreement with that 
found in the ASCA analysis of ICannon et alJ i 19991) and the 
more recent Chandra analysis of lMaAevItchetalJ J2000I) . but 
extends to larger radius (albeit with large uncertainties). 

A polytropic model, T °c , with the gas density de- 
scribed with the j6-model discussed above (Sect. 13.31 1. is a 
good description of the projected temperature profile. The best- 
fitting values for the central temperature and polytropic index 
are f = 8.09 + 0.23 keV and y = 1 . 15 ± 0.02, respectively; this 
model is also shown in Fig. [3] 

As a consequence of the high temperature, abundance mea- 
surements of this cluster are poorly-constrained, but are con- 
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Fig. 3. The projected radial temperature profile of A2218. Error bars 
are upper (lower) limits from fits with background normalisation at 
+ 10 (-10) per cent of nominal. The solid line is the best fitting poly- 
tropic model to the data. See text for details. 

sistent with being flat at the average value (~ 0.2) out to the 
detection limit. 



4.3. Temperature map 

We next produced a temper ature map using the a daptive bin- 
ning technique described in iBelsole et alJ feCXHal) . A full en- 
ergy band image with the sources excised was extracted from 
each camera, and the camera images coadded. Starting from 
the largest spatial scale (i.e., the entire image), the image was 
divided into cells. Each cell is divided by four until a certain 
threshold criterion is met (for this map, a total of 1000 counts in 
the full energy band prior to background subtraction). Spectra 
are then extracted from the source and background event lists 
with the resulting spatial binning, and background subtracted. 
The same absorbed MEKAL model was fitted to each spec- 
trum in the bands described above (this time with the abun- 
dance fixed to O.2Z ). The derived uncertainties are dependent 
on the temperature. In building the temperature map shown in 
Fig. [4] we only used cells with relative errors smaller than 50 
per cent. The temperature structur es in this map are ver y simi- 
lar to those in the Chandra map of lGovonietal]ll2004 . 



5. Mass analysis 

5.1. Deprojection of temperature profile and correction 
for PSF effects 

In order to have the best estimate of the mass profile it is 
necessary to take into account PSF and projection effects on 
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Fig. 4. Adaptively-binned temperature map of A2218. The contours 
are the same as those in Fig.Q 
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Fig. 5. The total mass profile of A2218, calculated from the X-ray 
data assuming hydrostatic equilibrium and spherical symmetry. The 
dashed line is the isothermal /3-model mass, assuming a constant tem- 
perature of kr = 6.6 keV (Sect. RITV The solid line is the best fitting 
isothermal sphere model (King approximation) for the total mass pro- 
file; also shown (dot-dash line) is the best fitting NFW mass model to 
the entire profile (see Table|2|for details). 



the temperature profile. Since correction for these effects us- 
ing simultane ous spectral fits under XSPEC tends to be unsta- 
ble JPointecouteau et al.Ll2004IPratt & Arnaudll2005l) we us e 
the method described in Pointecouteau, Arnaud & Prattl J2005t) . 
With the gas density represented by the best fitting yS-model 
(Sect. 13. 3> . we adopt a polytropic representation of the pro- 
jected temperature profile (Sect. 14.21 . and convolve this with 
a redistribution matrix which takes into account both PSF and 
projection effects. The projected profile was randomly sampled 
1000 times within the observed errors, the gas density profile 
being allowed to vary within the uncertainties of the best fit- 
ting yS-model. The final deprojected, PSF-corrected temperature 
profile and associated errors were then the mean and standard 
deviation of the 1000 corrected values at each point. 

Fitting a polytropic model to the deprojected, PSF- 
corrected temperature profile yields best-fitting values for the 
central temperature and polytropic index of to = 8.66 + 0.1 keV 
and y — 1.16 + 0.01, respectively. 

5.2. Total mass profile calculation 

Combining the deprojected, PSF-corrected density and temper- 
ature profiles (Sects. 1331 and I5TTV we can derive a total mass 
profile under the assumptions of hydrostatic equilibrium and 
spherical symmetry. Given the observed temperature structure, 
this is probably not a good assumption (certainly in the inner 
1 '.5 or so), but it is still interesting to see if X-ray and lensing 



mass estimates agree in the external regions. The mass pro- 
file was calculated usi ng the Monte Carlo method described in 
IPratt & Arnaud J2003I) : the resulting total mass profile is shown 
in Fig. |3 

For easier comparison with previous ana lyses, particu- 
larly results from ASC A JCannon et all ^999) and Chandra 
(Mac hacek et all Eo02). we also show the mass profile ob- 
tained when the temperature profile is assumed isothermal and 
the gas density is parameterised by a /?-model. For this we as- 
sume a temperature of kr = 6.6 keV (Sect. l4.lt and the best fit- 
ting yS-model in Sect. 13.31 The mass profile obtained using this 
method is slightly different from the one we have derived us- 
ing spatially resolved temperature information, being slightly 
lower in the inner regions, and higher in the outskirts. 

5.3. Mass profile modelling 

Since the mass profile is defined to very good precision, it 
is interesting to compare it with theoretical expectations. We 
thus tried fitting the mass profile with functional forms for the 
dark matter density, thes e being the King approximation to an 
isot hermal sphere (e.g., Sarazin, 1986) and the NFW profile 
frjavarro. Frenk & WhiteLll 997^ given by: 



M(r) = 4np c (z)6 c r 3 s m(r/r s ) (1) 



6 



G.W. Pratt et al.: An XMM-Newton observation of A2218 



Table 2. NFW and King isothermal sphere fits to the total mass 
profile of A2218. Errors are lcr. 



Parameter 


NFW Model 


King model 




Full 


Central point 
removed 




c 

r s (h~l kpc) 


4 96 +0 43 
345 + g 


c 70+0.60 

2S9 + _H 


q 74+0-52 
-0.50 

146+ 11 


rim (hjg kpc) 


1715 


1669 


1525 


M 2 m (M ) 


6.8 x 10 14 


6.3 x 10 14 


4.8 x 10 14 


x 2 /y 


14.2/8 


2.5/7 


4.8/8 



3 [ln(l + c) - c/(l + c)] 
m(x,NFW) = ln(l + x) - x{\ + x)' 1 (3) 
m{x, iso) = ln[x + (1 + x 2 ) 1/2 ] - x(l + x 2 y l/2 , (4) 

respectively, where we have followed 

lEttori. De Grandi & Molendil (|2002), in expressing the 
normalisation of both mass models through the concentration 
parameter, c. The King profile is clearly preferred in terms of 
X 2 , although the removal of the central mass point leads to 
an acceptable fit for the NFW mass model (see Table [2] for 
details). The best-fitting mass profile models are compared to 
the data in Fig.|5] 

6. Discussion 

6.1. Temperature, entropy and pressure structure 

It has been known for some time that the central regions of 
this cluste r show sign s of b e ing mo rphologi cally di sturbed 
in X-rav dMarkevitchl Il997l iNeumann & Bohringerl 1 19991 
iMachacek et all 120021) . The spatially r esolved temperature 
maps from this work and from Chandra JGovoni et all 120041) 
provide further evidence, if any were needed, that A2218 is 
undergoing (or has undergone) a merger. The centrally-peaked 
radial temperature profile has been resolved into a small region 
of significantly hotter gas in the core. The temperature varies 
by a factor of two in the central arcminute. The hottest parts 
of the cluster are only approximately spatially coincident with 
the weak radio halo dGovoni et alll2004T) . Comparing the HST 
galaxy distribution (Fig. Q and the temperature map (Fig. |4}, 
the hottest region is spatially coincident with the two main 
galaxies, and thus with the two main mass peaks found in lens- 
ing analyses, but is offset from the X-ray surface brightness 
maximum. 

We have also investigated the structure in entropy and pres- 
sure using masks created to resolve the observed fluctuations 
in (i) the temperature and surface brightness, and (ii) the en- 
tropy and pressure. Both masks have 16 regions, 8 regions cov- 
ering the central 1.'5, and 8 covering larger scales of 1' - 5'. 
The method of mask generation and the calculation of the de- 
projected pressure and entropy in each region is described in 
iHenrv. Finoguenov & BrieJ J2004I) . The fluctuations on these 
scales are found to be statistically significant and amount to 
(1 1 — 13) ±2 per cent in entropy and 14 + 3 per cent in pressure, 
where the two values given for the entropy show the slight de- 



pendence on the mask. Note however that the results obtained 
using the two different masks are consistent. The scatter in the 
entropy and pressure in the central arcminute is 17 + 3 per cent 
and 17 + 5 per cent, respectively. 

The temperature structure of a cluster has often been used 
to learn more about the merger kinematics in interacting clus- 
ters (e.g., Markevit ch et all fl999: Bels ole et alll2004albl) . This 
is very difficult in the present case, since, at least in X-ray, the 
merger geometry is not obvious. In particular, the X-ray sur- 
face brightness does not show more than one peak and is fairly 
symmetric at large radii. 

6.2. Entropy profile 

The entropy profile of the cluster shown in Figure [6] was 
calculated using the now-customary entropy definition, S = 

—2/3 

kTn e , where we have used the best-fitting jS-model for the 
gas density (Sect. 13 .31 and the deprojected, PSF-corrected 
temperature profile (Sect. 15 . II . Uncertainties due to gas den- 
sity profile modelling are negligible compared to the temper- 
ature measurement errors. Th e radius has been scaled u sing 
the ob served r2oo-T relation of lArnaud. Pointecouteau &" Pratt 
(2005J) 1 , and the entropy has been scaled using the relation 
found by Ponman et al. f d2003l) 2 . We also show an analytically- 
derived profile obtained using the best fitting analytic descrip- 
tion of the temperature and density data. 

It is clear that the entropy profile of A2218 is rather flat 
in the inner regions (r < 200 hz}. kpc, or 6 < 1'), where the 

— 1/3 

central entropy levels out to a value of S ~ 200 h J0 ' keV 
cm 2 at 20/i^J kpc. This flattening is driven by the large core 
in the gas density. Because of this, A221 8 was found to be a 
significant outlier in the sample of lLlovd-Davies et all J2000). 
In fact, the temperature map shows variations of a factor of two 
(from ~ 5 to ~ 10 keV) in the central arcminute. The complex 
entropy distribution in the central regions is smoothed out by 

the use of annular regions to produce t he profile. 

From an analysis of five clusters, IPratt & Arnaudl d2005l) 
found approximately 20 per cent dispersion in the scaled 
entropy profiles of their sample. Their mean scaled profile, 
obta ined in the ra dial range [0.05 - 0.5] r2oo (Eqn. 4 of 
Pratt & Arnaudl2005l) is also shown in Fig. [6] Excluding the in- 
ner three points, a power-law fit to the entropy profile of A2218 
yields a slope of a — 0.80 + 0.07, slightly shallower than the 
mean scaled profile, although consistent within the lcr uncer- 
tainties. 

6.3. Mass estimates 

6.3.1 . Mass profile: comparison to theoretical 
expectations 

It is becoming increasingly clear that the mass profiles 
of relatively relaxed, hot clusters are well described by 
the NFW profile (e.g., iDavid et all EoOll lAllen et alJ l200tt 

1 The relation is log r 200 = (2.81 + 0.02) + (0.57 ± 0.03) log T, giving 
r?oo = 1892 fci kpc for a global temperature of 6.6 keV (Sect l4.lV 

2 S cc h Z - 4l3 T°- 65 
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Fig. 6. The entropy profile of A2218. The radius has been scaled 
to r2oo using the rioo-T relation of lArnaud et alj $2005 ). The en- 
tropy is scaled using the empirical scaling of IPonman et alJ $2003, 
S oc h(zT 4/3 T 065 ). The line following the points is the entropy pro- 
file obtained from the best-fitting analytical models to the deprojected, 
PSF-corrected gas density and temperature profiles. The straight line 
is the mean scaled entropy profile found bv lPratt & Arnaud fc005l) . 
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Pointecouteau et al. 2004, 2005). Such clusters generally ex- 
hibit strongly peaked surface brightness profiles, and regular 
density and temperature structure. In cases such as these, the 
standard assumptions in the X-ray mass analysis (hydrostatic 
equilibrium and spherical symmetry) are more than likely war- 
ranted. 

An NFW model is not a good description of the total mass 
profile of A2218, particularly in the inner regions, as expected 
given the observed density and temperature structure. While 
exclusion of the central mass point improves the NFW fit con- 
siderably, the best fitting mass model to the entire profile re- 
mains the King model, a consequence more of the large core in 
the gas density than the radial temperature distribution. 

6.3.2. Comparison with lensing 

The considerable disagreement between the masses derived 
from X-ray and lensing anal yses of this cluster has generated 
some debate in the past (e.g.. lMiralda-Escude & Babullll995t 
lAllenUl998h . In Fig.[7]we show the projected mass profile ob- 
tained from these XMM-Newton data, converted to the SCDM 
(Qo = 10, Qa = 0.0, h = 0.5) cosmology used in previous 
analyses. The p rofile is compared to the masses obtained from 
strong lensing JLoeb & Mao[ll994t iKneib et all Il995t LMlenj 
1 19981) . weak lensing JSauires et all Il996f ISmail et all Il997h 
and combined lensing ( AbdelSalam et al., 1998; Smith etal., 
2004). 
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Fig. 7. The projected mass profile of A2218. Filled circles with er- 
ror bars are the XMM data points. The dashed line is the isothermal 
/?-model mass, assuming a constant temperature of 6.6 keV. Various 
other symbols indicate different lensing mass estimates. 



Despite having a well constrained radial temperature distri- 
bution, the X-ray mass is approximately two times less than the 
strong lensing mass at 85 h 5 ^ kpc. At larger radius, the agree- 
ment between mass estimates is better, although still not per- 
fect. The observed disagreement in the inner regions would be 
expected if either (i) there was a smaller subcluster in the line 
of sight, or (ii) the hydrostatic assumption is violated, possibil- 
ities we discuss in more detail below. 

6.3.3. Total mass estimate 

Depending on the adopted mass model, t he virial mass , M7.n0 , 
varies by about 30 per cent (see Table|2li. Arnaud et alJ (|2005) 
have recently derived an X-ray M—T relation derived from 
XMM-Newton observations of 10 clusters in the 2-9 keV 
range. Since the methods detailed in the present paper are very 
similar to their analysis, it forms an ideal comparison sam- 
ple. The tota l mass , M200, predicted by the M—T relation of 
lArnaud et alJ J2005I) . assuming a virial temperature of 6.6 keV 
(Sect. |47J3, is 8.5 x 10 14 h^M , about 25 per cent higher 
than the highest mass estimate from mass profile modelling 
(Table©. 

We can also compare the mass estimates from the differ- 
ent X-ray mass model s with the optical virial mass estimates 
of lGirardi &Mezzettil feOOll) . We note that the most distant 
galaxy from the cluster centre in Girardi & Mezzetti's anal- 
ysis is only at ~ 450 hz\ kpc, and so there is a fair amount 
of extrapolation needed in their mass estimate. These authors 



analyse 43 galaxies and find M V;0pt = (2.1 +0.5) x 10 15 M< 
for an optically-derived virial radius of 2.3 /i- ( J Mpc. Our best- 
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fitting mass models give 8.6/4.7xl0 14 /i TO ' M (NFW/King) for 
the same radius. This is a very significant mass discrepancy. If 
there is indeed a line of sight merger, this will artificially in- 
crease the optical galaxy velocity dispersion if only one com- 
ponent is assumed, leading to a higher optical mass estimate. 



7. Summary and conclusions 

We have undertaken a detailed X-ray analysis of the well- 
known lensing cluster A2218. In the inner regions (6 < V, 
or r < 200 fui kpc) the X-ray isophotes become elon- 
gated and oriented approimately SW-NE, the same direction 
as the line joining the two brightest cluster galaxies (which 
are also the mass peaks seen in lensing studies). There is a 
marked offset between the position of the cD galaxy and the 
X-ray brightness p eak, as was seen i n previous ROSAT and 
Chandra studies ( Markevitch, 1997; Neumann & Bohrinaer, 
1 19991 iMachaceketall 120021). In agre ement with previous 
Chandra results dMarkevitch etall EoOO). we find that the tem- 
perature profile decreases significantly from the centre. With 
the present observation, we find that this decrease continues 
to the virial radius of the cluster. From an adaptively-binned 
temperature map, we find similar temperat ure structure to that 
found in the Chandra observation jGovoni et al.Ll2004T) . There 
is a pronounced peak in the temperature in the inner 1', where 
the temperature rises from kT ~ 5 keV to kT ~ 10 keV. The 
peak is not correlated with any obvious surface brightness fea- 
ture, and is offset from the surface brightness maximum. The 
scatter in deprojected entropy and pressure in the central ar- 
cminute is typically ~ 17 ± 5 per cent. 

The mass profile (determined assuming HE and spherical 
symmetry) is best fitted with a King mass model, in contrast 
with the results from more relaxed clusters. An NFW model 
is a particularly bad fit in the central regions. The X-ray mass 
estimate does not agree with the mass estimated from strong 



lensing arcs at 80 /i 5( J , although the agreement improves when 
comparing X-ray and weak lensing masses at large radii. X-ray 
total mass estimates dependent on mass profile fitting can vary 
by 30 per cent depending on the mass model used. While the 
results must be viewed with caution, the total mass measure- 
ment is 20-60 per cent smaller than the mass expected from 
the M-T relation. X-ray and lensing mass measurements are 
significantly smaller than optical virial mass estimates by more 
than a factor of two. 

An important clue to the merger geometry in this system 
comes from optical observations . The distribution o f galaxy 
velocities in the optical study of Gira rdi et alJ l l 19971) strongly 
indicates substructure in the line of sight. Given this and the ob- 
served X-ray characteristics, it is likely that A2218 is either an 
ongoing merger between a main cluster and a smaller subclus- 
ter, or a single cluster in a more advanced merger state. In either 
case, the central regions, at least, of the cluster are unlikely to 
be in hydrostatic equilibrium. As a result of this, X-ray mass es- 
timates are highly uncertain. Discrepancies between X-ray and 
lensing mass estimates due to li ne of sig ht effec ts and/or merg- 
ers have als o been discussed inlAndersson & Madeiski (2004, 
A 1698) and Zhang et alJ J2004 CL0024+17), and in the case 



of A2218 itself, by iMachacek et alJ J2002I) . in the context of 
the Chandra observation. 

Perhaps the most significant result from our analysis is the 
calculation of a well-constrained radial temperature profile out 
to approximately the virial radius of the cluster. It is unfortunate 
in a way that A2218 is so obviously undergoing a merger event. 
In other circumstances it would have been an ideal benchmark 
object for mass measurements using the various methods in use 
today. 
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